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Abstract 

A general method to describe a second-order phase transition is discussed. 
It starts from the energy level statistics and uses of finite-size scaling. It 
is applied to the metal-insulator transition (MIT) in the Anderson model of 
localization, evaluating the cumulative level-spacing distribution as well as 
the Dyson-Metha statistics. The critical disorder W c = 16.5 and the critical 
exponent v = 1.34 are computed. 
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The Anderson model is well known to ad- point W c , 
equately describe some characteristic aspects 

u(w) = c\w-w c r. (i) 

of disordered systems. In three dimensions 

(3D) it exhibits a metal-insulator transition Here W is the critical parameter and ( a 

(MIT). In spite of the simplicity of the An- proportionality constant. In the usual TMM 

derson Hamiltonian, the critical exponent re- the correlation length £m of a finite system 

mains a controversial subject theoreticallyBH of size M is calculated and extrapolated to 

as well as numerically.1! Recently previous an infinite system size by a finite-size scaling 

contradictions concerning the numerical re- procedure. In the following a characteristic 

suits could be removed! and a universal value quantity olm is used from which the critical 

for the exponent established for several nu- parameters are derived. There is no need to 

merical data sets. But all these numerical compute £m and one has a large choice for 

studies were done in the framework of the the definition of 

transfer matrix method (TMM).i Here we The present paper extends previous 
discuss another approach, no longer based on studies^ of the statistical properties of the 
the TMM data, to calculate critical parame- energy spectrum of the Anderson Hamilto- 
ters like the critical disorder W c and the crit- nian on both sides of the MIT. In a similar 
ical exponent v and confirm the previous re- way the influence of the magnetic field on the 
suits. The approach is fairly general and can spectrum with! or without!! disorder has been 
be easily applied to other phase transitions, considered. After unfolding the spectrum0 
Only for clarity of our arguments we formu- its fluctuations are characterized by means of 
late the subsequent derivation in terms of the the spacing distribution P(s) and the Dyson- 
Anderson model of localization and its MIT. Metha statistics A3. The first entity mea- 
The MIT is expected to be a second or- sures the level repulsion between consecutive 
der phase transition which can be character- levels and the second one reflects the spectral 
ized by the correlation length ^ of the order rigidity. In terms of these quantities the MIT, 
parameter fluctuations. This length diverges in a system of finite size M, is a continuous 
with the critical exponent v at the critical transition from the Gaussian orthogonal en- 



semble (GOE) to the Poisson ensemble (PE). for a finite size. As a(M, W) describes 
These two limiting regimes, GOE and PE, a system which exhibits a second order 
were derived by Alt'shuler et al. using dia- phase transition with characteristic length 
gram technique for conduction electrons mov- £oo, Eq.(l) suggests generally, following the 
ing in a random impurity potential and tested work of Brezinjll a finite-size scaling be- 
numerically for a small system (M = 5).0 haviour for a(M, W), namely a(M, W)/a = 
The values of P(s) and A 3 are functions of g^M/^iW)) or a(M,W) = f (M / ^(W)) . 
the size M except at the critical point. In But as a(M, W) is analytical for finite M, / 
the thermodynamic limit this critical point can be expanded around the critical point as 
behaviour is expected to separate two differ- f(x) x =° a + bx x l v + cx 2 ^ u + . . ., yielding 
ent regimes: GOE for the metallic side and a(M, W) ~ a(M, W c ) + C\W — W c \M l ^ u 
PE for the insulating side. Similar results 

(3) 

have already been mentioned by Evangelou et 

aM and studied in more detail by Shklovskii for an Y quantity a which fulfills the condition 
et a/I but considering only P(s). We con- ( 2 ) with a singularity at W c for M -»• oo and 
sider a yet undefined characteristic quantity which obe y s the scaling relation a(M, W) = 
a{M, W) which shall be chosen in such a way f(M/U(W)). We note that one has to con- 
that: sider a finite-size scaling0 procedure so that 

q, yjy < the size of the system is kept finite, suitable 

yY = W c (2) f° r numerical work, and not usual scaling^ 

- yyy > for which an infinite system has to be consid- 

where a Me, otin and occr ar e different con- ere d. 

stants for the metallic side, the insulating Before defining a(M, W) we recall how 

side and at the critical point, respectively. p ( s ) and A 3 were computed! The energy 

With respect to the W dependence a (M,W) spectra were obtained from the Anderson 

is singular for M = oo, but analytical Hamiltonian H = £i e f |z)(z| + E^j V\i)(j\, 

otherwise because all the singularities con- where the sites % are distributed on a sim- 

nected with the phase transition are rounded P le cubic lattice, interactions with only the 



nearest neighbours are considered and V=l 
defines the energy scale. The site energy e.; 
is described by a stochastic variable given by 
a box distribution of width W which repre- 
sents the disorder and is our critical param- 
eter. The secular matrix of the Hamiltonian 
was diagonalized in a straightforward way by 
means of the Lanczos algorithm for systems 
of size M 3 from 13 3 to 21 3 and disorder W 
from 2 to 40. The number of different re- 
alizations of ti was chosen so that ~ 2 • 10 5 
eigenvalues were obtained for every pair of 
parameters (M, W). Thus the full spectrum 
was computed 25 to 90 times. For the subse- 
quent investigations only half of the spectrum 
around the band center is considered.0 

To analyse the scaling of the statistical 
properties of the eigenvalue spectrum, start- 
ing from P(s) an evident choice for a, which 
satisfies the expected finite-size behaviour, 
would be a = J °° P(s)ds. But because P(s) 
is normalized this makes no sense. A solu- 
tion is provided by the appearance of a fixed 
point at So — 2 when one plots P(s) for dif- 
ferent disorders W and sizes M.i'B So one 
can choose ap( s )(M,W) = P(s)ds, with 
so — 2 as already quoted.! An equivalent 
choiceB would be a P{s) (M,W) = f°° P(s)ds. 



While these choices are formally possible it is 
much better from a numerical point of view 
to consider the cumulative level-spacing dis- 
tribution J(s) = Jq P(s')ds' as plotted in 
Fig.[TJ instead of P(s) because the most signif- 
icant changes which occur in P(s) for small 
s are emphasized and because I(s) is much 
smoother. Using I(s) we choose 

1 r s o 

a x (M,W) = — / I(s)ds . (4) 
s Jo 

The results in FigJ^ clearly distinguish three 
different regions, in which a% increases, re- 
mains constant, or decreases with increas- 
ing size M, respectively. This is consis- 
tent with the properties (2) of a(M,W) re- 
quired in order to apply the method. We 
can already determine the critical disorder 
from Fig.||: ai\ (M, W) is size independent for 
W c = 16.25 ±0.25. 

In the next step we check whether 
«i(M, W) can be expressed by a scaling func- 
tion f{M/^ s00 (W)). To perform this scaling 
procedure with the a,\ data only, the range 
of ai values for various M at any given dis- 
order W\ must overlap the range of a\ values 
for various M for at least one different dis- 
order W-2- That would require to compute 
oti(M,W) for more W values or larger M. 



This is beyond our present computational 
means. Because of this restriction we use the 
values of ^ from previous calculations done 
in the framework of the TMM.0 Of course, 
this does not contradict the statement in the 
introduction that the method is no longer 
based on the TMM; only the present eval- 
uation is. The TMM data are not necessary 
for the calculation of W c and v below but just 
for the demonstration of the scaling. The re- 
sult in Fig|| confirms the scaling hypothesis 
for «i(M, W). The upper branch of the curve 
represents W > W c and the lower W < W c . 
The fluctuations in the scaling curve can be 
explained by the facts that one uses a cor- 
relation length obtained from a different 
method, that some values of ^ in the crit- 
ical region had to be interpolated, and that 
the range of the a\ values is very small. The 
biggest fluctuations are due to the smallest 
M which indicates that M = 13 may be not 
large enough or that one has to average over 
more realizations. 

In the final step we determine the crit- 
ical exponent v. The tricky point is that 
the formula (3) is valid only in the vicin- 
ity of the critical point but on the other 
hand one knowsifl that close to the transition 



the numerical inaccuracies are largest. The 
task is therefore to find an adequate range 
for \W — W c \ which satisfies these two con- 
straints. This can be controlled by a \ 2 ~ 
fittingS which measures the agreement be- 
tween the data and the "model" (3). To de- 
termine the "best" range we minimize x 2 /N ', 
where N is the number of data considered, as 
a function of the range \W — W c \ and obtain 
v = 1.37±0.15 for \W — W c \ ~ 1 in very good 
agreement with our previous resultS v ~ 1.34 
for the box distribution. Using «p( s ) but with 
even fewer data around the critical point and 
for smaller samples, Shklovskii et a/.i esti- 
mated v ~ 1.5 which is consistent with our 
results. 

Considering as above I{s), only a part of 
the information contained in the spectrum 
is evaluated, namely the correlation between 
consecutive levels. By the A 3 -statistics (see 
Figl4|) the correlation between nonconsecu- 
tive levels is taken into account, too. Now 
we choose 

1 r 30 

a 2 (M,W) = - A 3 (L)dL (5) 
3U jo 

The results in Fig.[5] show that the critical 
disorder can be located at W c = 16.5 ± 0.1. 
The quality of these data is better than that 



of «i(M, W): the fixed point is well defined crepancy between the numerical predictions 

and the different M dependencies described for v and the experimental values^ between 

above for a\ can be seen more accurately for 0.5 and 1, raised the question whether the 

CK2- The scaling assumption for a2(M, W) is Anderson model is relevant to describe dis- 

tested in the same way as that for a\ and ordered systems. But recently experimental 

confirmed in Fig|| For the critical exponent, measurements of v in uncompensated Si:P 

using Eq.(3) and the x 2 fitting, one obtains samples have yielded v = 1.3.0 So while 

v = 1.34 ±0.10, identical to the value previ- it cannot yet be claimed that the problem 

ously found,@ with \W — W c \ ~ 1. As sup- of the critical exponent in disordered sys- 

posed above these results are more accurate terns is solved it seems that results obtained 

than the ones provided by P(s) or J(s). in very different ways converge towards the 

In conclusion, the here employed method same value v ~ 1.35. 
to study the MIT on the basis of the statis- It is surprising that the same value has 
tical properties of the energy spectrum pro- been found applying the TMM to 3D sys- 
vides very good results qualitatively as well terns with magnetic field.B The analysis of 
as quantitatively. The obtained critical expo- the statistical properties of the energy spec- 
nent v = 1.34 ± 0.10 confirms the previous^ trum can be a promising way to understand 
TMM result. Applying the same TMM but a this. Until now it was claimed that one could 
slightly different way to compute the critical expect different critical exponents for differ- 
exponent from the raw data, MacKinnonll! ent universality classes e.g. the orthogonal 
finds v ~ 1.53. While that result is not in (GOE) and the unitary (GUE) classes. How- 
contradiction with our uncertainty, we be- ever, new results have shownUll this type of 
lieve that the difference is significant. The- symmetry only for small disorder but not at 
oretically the derivation of the critical ex- the critical point where one finds an inter- 
ponent remains an open subject although mediate distribution e.g., between GOE and 
recently Hikami,^ using string theory, pro- PE in the orthogonal case. The same critical 
posed the possibility to obtain v = 4/3 which exponent should mean the same universality 
would be supported by our results. The dis- class at the critical point. Thus the observed 



agreement of the critical exponents with and 
without magnetic field may be interpreted as 
a signature of a new universality class differ- 
ent from PE, GOE, and GUE.i 

Finally we stress that although the 
method was employed for the MIT in the An- 
derson model of localization, the approach is 
very general and easy to apply. It would be 
very interesting to study other second order 
phase transitions from this point of view not 
necessarily considering the eigenvalues of the 
Hamiltonian, but for example the eigenvalues 
of the transfer matrix in Ising models. 
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FIGURES 

FIG. 1. Cumulative level-spacing distribu- 
tion I(s) for a system size 21 x 21 x 21 and the in- 
dicated values of the disorder W. The solid lines 
represent I(s) for the GOE and the PE cases. 

FIG. 2. Scaling variable a\{M,W) derived 
from I(s), according to Eq.(4). For the GOE 
and the PE the limiting values &Me — 0.507 and 
&in — 0.569, respectively, can be derived. 

FIG. 3. Verification of the scaling assump- 
tion for a%(M, W). The correlation length ^ is 
taken from TMM dataH 

FIG. 4. A3-statistics for a system size 
21x21x21 and the values of the disorder W 
indicated on the right hand side. 

FIG. 5. Scaling variable a<i{M, W) derived 
from A^(L) according to Eq.(5). For the GOE 
and the PE the limiting values &Me — 0.264 and 
5/n — 0.804, respectively, can be derived. 

FIG. 6. Verification of the scaling assump- 
tion for ct2(M, W). The correlation length ^ is 
taken from TMM data!! 
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